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ACTIVI A: Calculation of isotope production cross-sections and yields 

J. J. Back*, Y. A. Ramachers 

Department of Physics, University of Warwick, Coventry, CV4 7AL, United Kingdom 



Abstract 

We present a C ++ computer package, ACTIVIA, that can calculate target-product cross-sections and the production 
and decay yields of isotopes from cosmic ray activation using data tables and semi-empirical formulae. We describe 
the structure and user interface of the computer code as well as provide comparisons between the calculations and 
experimental results. We also outline suggestions on how the code can be improved and extended for other applications. 
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1. Introduction 

There are many observations made in astro- 
physics, astronomy and particle physics that require 
knowledge of the cross-sections of nuclear processes 
at high energies. The physics of cosmic rays is one 
example. Understanding how they interact with 
matter allows us to study the composition of the in- 
terstellar medium, the Galaxy and the solar system. 
Estimating the cosmic ray background radiation 
in spacecraft and isotope production in accelerator 
facilities are other important applications. For the 
case of low background experiments, it is necessary 
to know which isotopes are formed by cosmic ray 
activation inside target materials before they are 
shielded underground. Any long-lived isotopes must 
be accounted for in the signal data analysis, and 
this requires information about cross-sections and 
production rates. There are measurements available 
for use in such studies, but normally only for a re- 
stricted energy range which does not cover the total 
energy spectrum of the input beam. Some impor- 
tant reactions for low background experiments, such 
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as Co produced from natural tellurium by cosmic 
ray activation, have essentially no experimental 
data available. Accordingly, a set of semi-empirical 
formulae have been developed by Silberberg and 
Tsao [112] to estimate the cross-sections of various 
nuclear processes, such as spallation and fission, 
which can be used to calculate the production yields 
of nuclear isotopes from cosmic ray activation. The 
parameters in the formulae, assumed to be the same 
for proton and neutron beams, have been tuned to 
best fit the available experimental results. 

In this paper we describe a new C ++ package, AC- 
TIVIA, that uses a combination of semi-empirical 
formulae and tables of data based on experimen- 
tal results to calculate the cross-section, production 
rates and yields of radioactive isotopes from cosmic 
ray activation. Section [3] describes the Silberberg- 
Tsao semi-empirical formulae used for the cross- 
section calculations while Section [3] explains how 
cross-sections can be evaluated if tables of data (in 
ASCII text format) are available. Section[4]provides 
details on how the radioactive decay yields are cal- 
culated, Section[5]explains the basic structure of the 
code, while Section [5] shows comparisons between 
the calculations and experimental data. Finally, in 
Section[7]we describe possible modifications and ex- 
tensions that could be made to the code if other in- 
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put beam spectra or nuclear physics processes are 
required. 

2. Semi-empirical formulae 

We use the extensive set of semi-empirical for- 
mulae from Silberberg and Tsao |l|2j to calculate 
cross-sections for isotope production when a proton 
or neutron beam (e.g. cosmic rays) hits a given tar- 
get. The general form of the cross-section equation 
for a product (Z, A) from a target (Z t , A t ) is 

a = a Q f(A)f(E)e- pAA x (1) 
exp (-R\Z -SA + TA 2 + UA 3 \ V ) Slrfi, 

where E is the beam energy (in units of MeV), A A is 
the difference between the target and product mass 
numbers (A t — A), while f(A) and f(E) are correc- 
tion factors usually applied for heavy targets (A t > 
30) and when A A is large (> 10). The normalisation 
factor ctq is the cross-section (in millibarns) at the 
energy Eq (~ 3 GeV), above which cross-sections are 
assumed to be independent of energy. The first ex- 
ponential term describes the reduction of the cross- 
section as AA increases, where P is a function of tar- 
get mass At and E. The second exponential term de- 
scribes the distribution of cross-sections for several 
isotopes for a given element of atomic number Z; its 
parameters (R, S, T, U and v) are, in general, func- 
tions of target and product Z and A values as well 
as energy. The parameter R represents the width of 
the cross-section distribution amongst the isotopes, 
while S describes the peaks for various mass num- 
bers A. The parameter T describes the increase in 
the number of neutrons as A increases, while U is 
equal to3.0xl0 -7 and represents a small correction 
to improve the overall agreement with available ex- 
perimental data. The power exponent v varies from 
1.3 to 2.0 depending on the target and product iso- 
topes. The parameter fi represents a nuclear struc- 
ture factor, rj is a nuclear pairing factor (for even 
and odd Z and neutron (JV) numbers), while £ is 
an enhancement factor for light evaporation prod- 
ucts; these parameters are of order unity. We assume 
that Eq.[T]is applicable to both proton and neutron 
beams, with identical formulae and parameters. 

Figure [1] shows the various regions for different 
hadronic processes which can be calculated using 
Eq. [1] and which are implemented in the computer 
code. See Silberberg and Tsao |H2j for detailed in- 
formation about the parameterisations of the var- 
ious factors in Eq. [TJ The code also calculates the 
cross-sections for tritium production using the for- 
mulae from Konobeyev and Korovin [3] . 
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Fig. 1. Domains of applicability of different parameters of 
Eq. U] as a function of target mass At and product mass 
A, which are implemented in the code. P denotes the re- 
actions that are mainly peripheral; SL, SI and SH denote 
the spallation of light, intermediate and heavy targets, re- 
spectively. The fission region (including contributions from 
spallation reactions) is denoted by F, the symbol B repre- 
sents the breakup region in which light nuclei are evaporated, 
while BF represents the combination of fission and breakup 
reactions. The evaporation of light nuclei occurs in region 
e, with Se combining evaporation and spallation processes. 
These are the same regions as those in Ref. pp. Finally, the 
dotted line represents the tritium production region 

3. Data tables 

There are extensive sets of nuclear data (e.g. the 
EXFOR database [4]) that can be used to provide 
cross-sections as a function of energy for target- 
product isotope pairs. If present, these (ASCII text) 
tables are used instead of the semi-empirical for- 
mulae. A linear interpolation is performed between 
data points to provide the cross-section (in mb) at 
a given beam energy (in MeV). It is important that 
the data be in the correct units, have enough energy 
bins to be useful and have a contiguous, increasing 
energy range. Equation Q] is used to calculate the 
cross-section if any data value is below a selectable 
lower limit (recommended to be set at 0.001 mb). 
If two or more separate energy regions are needed, 
then the data from these must be combined into the 
same table. For the case of large energy differences 
in between neighbouring regions, it is sensible to set 
the cross-sections to zero just above and below the 
range limits in the table so that the semi-empirical 
formulae are used instead. Otherwise, the linear in- 
terpolation would only give a weighted average of 
the cross-sections between the data limits of neigh- 
bouring energy regions, which may not be reliable 
enough. The code can use MENDL-2P tables [5] for 
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intermediate energy (up to 200 MeV) cross-sections. 



F/ xp = Yi (1 - e - x ^) , 



(4) 



4. Radioactive decay yields 

Consider a target comprised of several isotopes 
with the same atomic number Z t but having differ- 
ent mass numbers An each with a relative abun- 
dance fraction fa such that Ylifi ~ 1- This tar- 
get is exposed to a beam with an energy spectrum 
d<f)/dE (cm -2 s -1 MeV -1 ), such as cosmic ray neu- 
trons shown in Fig. [2j which is based on the param- 
eterisation from Armstrong \<$\ and Gehrels [7] . The 




Energy (MeV) 



Fig. 2. Energy spectrum for cosmic ray neutrons at sea 
level based on the parameterisation from Armstrong [6] and 
Gehrels \7\. 

production rate (kg -1 day -1 ) of an isotope j will be 

where C is a normalisation factor and cr^ (E) is the 
cross-section (mb) for the target-product isotope 
pair at energy E (MeV). The integral is over the 
required energy range of the input beam spectrum, 
while the sum is over all of the target isotopes. Here, 

C = N A x 24 x 36 00 x 10 - 24 , (3) 

where Na is Avogadro's number (6.022 xl0 23 mol -1 ), 
the next two numbers convert the production rate 
from per second to per day, while the last factor 
accounts for the different units of area used for the 
flux (cm 2 ) and cross-section (mb) quantities, as well 
as the unit of mass used for the atomic number An 
(converting g to kg) . If the target is exposed to the 
(continuous) beam for a time i exp (days), then the 
yield rate (kg -1 day -1 ) of the radioactive product 
isotope j with a half-life of t 1 / 2 (days) is given by 



where Xj = ln2/f 1 y 2 and assuming no chained de- 
cays. If the target is removed from the beam and the 
product isotope is allowed to decay for a time tdec 
(days), then the final product yield rate is 

ydcc = y^c X p e -A jtdoc = y. ^ _ e -A 3 i oxp ) e -A 3 W_( 5 ) 

Here, we only consider product isotopes with half- 
lives greater than one day. However, there are many 
short-lived isotopes which can decay to these, and 
they are included in the cross-section and yield cal- 
culations as "side-branch nuclei". For example, con- 
sider the decay chain to the long-lived isotope g 4 C 
(t 1/2 = 5730years): | 4 Be -> ^ 4 B +e - + D e (t 1/2 = 
4 ms) and £ 4 B -» 14 C +e - + v e (t 1/2 = 14ms). This 
means that both 4 4 Be and g 4 B are considered side- 
branches of | 4 C, and their cross-sections and yields 
will contribute to those for | 4 C. 

5. Code structure 

The code, written in C ++ to take advantage of 
object-oriented design principles, is split into two 
main sections; the first part performs the calculation 
of the cross-sections (Eq. [I} and production rates 
(Eq. [2]), while the second part calculates the final 
yields of the radioactive isotope products (Eq. [5]) . 
The basic structure of the code, designed to work in 
the Unix environment, is described below. Details 
on how to build and run the code are found in the 
README file of the computer package [8] . 

5.1. Input 

There are several inputs that are necessary for 
the cross-section and yield calculations. They are 
the target isotopes, a list of all allowed product iso- 
topes, the input beam spectrum, as well as the beam 
exposure and target cooling times. These are spec- 
ified by the ActAbsInput abstract interface, which 
is a base class that contains several pure virtual 
functions that must be implemented by any derived 
classes, such as Actlnput. Figure [3] shows a simpli- 
fied version of the overall program flow. 

The class ActTarget defines the target element, 
together with its isotopes and their relative abun- 
dance fractions /, (such that J2i fi = !)• If a target 
is comprised of many elements (each with their own 
isotopes), then the code must be run as many times 
as is necessary to get the results for each target el- 
ement. The results must then be combined at the 
end by the user. 
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Fig. 3. Diagram showing a simplified version of the program flow and C + + class structure. The input specifies the beam 
spectrum, the list of product and target isotopes, the decay times, as well as algorithms for the cross-section and yield 
calculations. These algorithms use the various inputs with the required cross-section and decay yield models to generate the 
output files for further analysis. 



The list of allowed products is denned in the Act- 
ProdNuclideList class, which is built from an ASCII 
text file ("decayData.dat") containing a complete 
list of radioactive isotopes found in the periodic ta- 
ble with half-lives equal to or greater than one day. 
All target and product isotopes are represented by 
ActNuclide objects, which store the atomic (Z) and 
mass (A) numbers as well as the half-life ii/ 2 - Prod- 
uct nuclide objects also contain the Z and A values 
of any side-branch nuclei (defined at the end of Sec- 
tion[3J . Each nuclide is created only once by the fac- 
tory class ActNuclideFactory, independent of how 
many times it is used in the code. 

The input beam spectrum is defined through the 
abstract class ActBeamSpectrum, in which the flux 
(d(f)/dE in units of cm~ 2 s" 1 MeV -1 ) must be im- 
plemented by any derived class. At present, the only 
available input beam is the cosmic ray spectrum 
shown in Fig. [2] which is assumed to be valid for neu- 
trons and protons. The ActBeamSpectrum class also 
stores the energy range and bin width A.E, which 
are used to calculate the integral in Eq. [3J as well as 
the atomic number Zbcam and mass number Ab ea m 
of the beam particles, both equal to unity for cosmic 
ray protons. 

The final required input variables are the time t exp 
that the target is exposed to the beam and the time 
elapsed after the target is removed from the beam 
(decay or cooling time idec)> both in units of days. 
These are stored within the simple ActTime class. 

The set-up of the calculations is managed by the 
Actlnput class, as mentioned earlier, which prompts 



the user to provide the required inputs using the 
standard C ++ input streamer "cin" . Examples are 
provided in the README file of the package [8] . 

The algorithms for calculating the cross-sections 
and yields must also be selected in the input class. 
At present, the cross-sections are calculated using 
Eq. [T] (with any available data tables as described in 
Section[3]) by the ActSTXSecAlgorithm class, while 
the yields are evaluated using Eq. [5] by the Act- 
SimpleDecayAlgorithm class. We now describe how 
these algorithms work. 



5.2. Cross-sections 

The calculation of the cross-sections between all 
target and product isotope pairs is controlled by the 
ActlsotopeProduction class, which retrieves infor- 
mation about the target, the list of product isotopes, 
and the beam energy spectrum from the input class. 

The cross-sections for each target isotope (Z t , A t ) 
are calculated using the ActProdXSecData utility 
class. This loops over all available product isotopes 
(Z, A) and calculates the cross-section for each en- 
ergy bin in the input beam spectrum. The values of 
Z and A must satisfy the requirements Z < Z t + 
Zbeam and A < A t + ^4bcam — 1, which ensure that 
there are at least two isotope products in each re- 
action, i.e. the cross-sections for any + p —>z+i 
Y + 7 reactions are not calculated. Any side-branch 
nuclei are included in the cross-section calculations. 
The energy of the reaction must be larger than a 
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threshold energy, typically of order 10 MeV, defined 
as the mass excess from the target-product reaction, 
and is calculated using the semi-empirical formula 
from Seeger [9]. Otherwise, the cross-section is set 
to zero and the next energy bin is considered. This 
requirement avoids calculating the cross-section of 
resonances for low-energy (E < 10 MeV) reactions, 
which is beyond the scope of this work. 

The required cross-section is retrieved from the 
ActSTXSecAlgorithm class. This first checks if there 
are any cross-section data for the target-product iso- 
tope pair at the given energy (above threshold). If 
no such data exists, or the cross-section is too small 
(< 0.001 mb, for example), then Eq.[T]is used to cal- 
culate the cross-section. The algorithm selects which 
calculation region in Fig.[T]is appropriate, i.e. what 
parameters are needed for Eq. [TJ and calls the class 
that implements the specific Silberberg-Tsao model, 
such as spallation, fission or peripheral reactions. 
The total cross-section is then just the sum of a{E) 
over all energy bins. The calculations do not distin- 
guish between ground or metastable states for the 
product nuclei. 



5.3. Yields 

The production rates of the isotopes (Eq. [2]) are 
calculated by ActlsotopeProduction when the cross- 
sections are evaluated. Following Eq. [2J the produc- 
tion rate at the energy E^ is given by 



yj (E k ) = Cj2j^%<nj(E k )&E, 



(6) 



where AE is the energy bin width from the input 
beam spectrum. The total production rate is given 
by the sum of yj(Ek) over all energy bins 



Y, 



A 

A t , 



J2^ aij (E k )AE 



(7) 



where the sum over k specifies the required energy 
bin range. The ActSimpleDecayAlgorithm class im- 
plements Eq. [5] to find the final yield Y^ ec for each 
product isotope j from the given exposure and decay 
times. This simple calculation ignores any chained 
decays, in which isotopes can decay to several other 
long-lived nuclei, or when isotopes can be made from 
multiple decays. However, short-lived isotopes are 
included in the final product yields through the use 
of the side-branch nuclei, as mentioned in Scction[3J 



5.4. Output 

The results of the cross-section and yield calcula- 
tions are stored in output files using the abstract in- 
terface ActAbsOutput, which contains several pure 
virtual functions that must be implemented by any 
derived class, such as ActStreamOutput for ASCII 
text output or ActROOTOutput for ROOT [T^ out- 
put. These classes must be able to deal with three 
different formats: lines of text (character strings or 
ActString objects), tables of numbers (ActOutput- 
Table objects) and graphs (ActAbsGraph objects). 
The cross-section and yield algorithms are not con- 
cerned with specific formats of output files; they only 
create tables and graphs which are then passed onto 
the output class, which decides how to write out the 
information in whatever format is required for fur- 
ther analysis. 

The ActOutput Tables store rows of numbers, 
where each column represents a specific variable 
with its own unique label. Graphs are made up of 
a series of points, which are represented by Act- 
GraphPoint objects that can store several y values 
for a given x co-ordinate, each with their own la- 
bels and units. For example, ActXSecGraph has x 
equal to the energy E (MeV) , y\ equal to the cross- 
section value (mb) and y 2 equal to the production 
rate Yj (kg -1 day -1 ). ActDecay Graphs contain 
points storing the elapsed time (in days) and the 
corresponding isotope yield: [x, y) = (t exp , V, exp ) 

or (tcxp + idcc, yf cc )- 

Each output class must be able to deal with two 
output files; the first one contains details from the 
cross-section and production rate calculations, while 
the second one contains the yields from the isotope 
decay algorithm. 

The total cross-section (mb) and production rate 
(kg -1 day -1 ) are stored in the cross-section output 
file. Also stored is the cross-section and production 
rate as a function of input beam energy (MeV), but 
only if a flag controlling the level of output detail 
is set appropriately. Again, the output class must 
determine what to do when the level of detail is 
changed; the cross-section and yield calculations are 
not affected by this. Finally, a summary table of the 
yields of each product isotope at the start and end of 
the decay period (t cxp and t cxp + idee) , weighted over 
all target isotopes, is written to the decay output 
file. 
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Fig. 4. Cross-section as a function of energy for (a) 60 Co and (b) 128 1 produced by protons in natural tellurium. Data sets in 
(a) are from Norman [12] and CERN [13] . Data points in (b) are from Acerbi [14j . Kormali [15] and Scholten [16] . 




Fig. 5. Cross-section as a function of energy for (a) 68 Ge and (b) 6t) Co produced by protons in natural germanium. Data points 
are from Norman | 12| . Horiguchi 1171 and Aleksandrov |18| . 



6. Comparisons between measurements and 
calculations 

The calculations of the cross-sections have been 
benchmarked using experimental data from the EX- 
FOR database [3] for natural tellurium, germanium 
and copper targets. The results of the calculations 
agree in general with those shown in various reports 
from the IDEA collaboration [11] . Figures H] to [6] 



show the calculated cross-section graphs of several 
isotopes along with the data points from experimen- 
tal results [3). The agreement between calculations 
and experimental results is not uniform; some re- 
sults agree very well (<50% relative difference), oth- 
ers agree within a factor of two, but there are also 
discrepancies up to a factor of 10 or more. It is clear 
that some of the MENDL-2P [5] cross-sections for 
energies between ~10 MeV to 200 MeV do not agree 
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Fig. 6. Cross-section as a function of energy for (a) 65 Zn and (b) 60 Co produced by protons hitting a natural copper target. 
The calculation of the 60 Co cross-section for energies in the range 40-200 MeV uses the MENDL-2P data tables, which 
underestimate the values obtained by various experiments. Data points shown in (a) are from Takacs 1191 . Michel'97 |20| . 
Szelecsenyi [21] , Kuhnhcnn 22 , Schiekel [23) . Cumming |24) and Michel'95 |25| . Data points shown in (b) are from Ozafran [26] ■ 
Takacs [S7|, Greenwood [28], Titarenko [29], Sisterson [30], Kim [31], Griitter [32], Orth [33], Aleksandrov [34], Michel'86 [35], 
Heydegger [36], Cline [37] and Mills [38]. 



very well with experimental results, such as those 
shown in Fig. [SJ Obviously, the agreement will im- 
prove if the MENDL-2P data tables are updated ac- 
cordingly. In fact, the user is free to use other exper- 
imental data for the cross-section tables instead of 
those from the MENDL-2P database. Instructions 
on creating self-made data tables are provided in the 
README file of the package [§]. 

Table 1 

Induced activities y dcc (^Bqkg" 1 ) of Te02 for three nu- 
clides assuming four months of exposure to cosmic rays and 
two years of decay in a shielded underground environment. 
There are two calculations from the IDEA Collabora- 
tion IDEA a uses the cosmic-ray flux distribution </> from 
Fig. l2l while IDEA b uses the distribution from Ziegler 1391 . 





60 Co 


125 Sb 


124 Sb 


CUORE [40] 


0.20 


15 


0.05 


IDEA a rn] 


0.02 


6±6 


0.01 ±0.03 


idea 13 rrn 


0.02 


22 ± 15 


0.05 ±0.06 


ACTIVIA 


0.03 


6 


0.02 



Comparisons have been made between the in- 
duced activity calculated by ACTIVIA with results 
shown in reports from the IDEA collaboration [11] 
for several targets and isotope products. Table Q] 
shows an example of the calculated induced activi- 
ties for Te02 exposed to cosmic rays at sea level for 



four months and shielded underground for a further 
two years. Note that the ACTIVIA result for 60 Co, 
which predicts a production rate Yj = 0.8^Bqkg _1 
(YIELDX [41] predicts OJ^Bqkg- 1 ), is a factor of 
10 lower than the estimate obtained by the CUORE 
collaboration [40] (using a modified version of 
COSMO [12]) owing to the improved calculation of 
the cross-section (see Fig. [4j. Tables [2] and [3] show 
the induced rates of isotope production by cosmic 
rays at sea level for natural and enriched germa- 
nium, respectively. One can see that there are dif- 
ferences between the various estimates presented in 
the tables. For example, the ACTIVIA results using 
the MENDL-2P data tables shown in Table [2] are 
generally quite close to those from Avignone [43] . 
but are approximately a factor of two less than those 
from HMS- ALICE @1] and YIELDX gT] estimates 
and those from Klapdor-Kleingrothaus [45] . Similar 
discrepancies between measurements and calcula- 
tions are also found for the enriched germanium 
yield estimates shown in Table [3] Table [4] shows 
the isotope production yields for natural copper. 
Two of the ACTIVIA results agree with those from 
Baudis and Schnee [49], but all of the calculated 
yields are less than those estimated by the IDEA 
group, whose calculations included the MENDL low 
energy (E < 10 MeV) neutron data [50]. Most of 
the cobalt results from Formaggio [IB] in Tables 
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Table 2 

Induced rates of isotope production Yj (kg -1 day -1 ) in natural germanium at sea level. The values from Avignone I43| are 
separated into Monte Carlo simulation (MC) and experimental (exp) results. 





68 Ge 


60 Co 


65 Zn 


58 Co 


57 Co 


56 Co 


54 Mn 


63 Ni 


55 Fe 


HMS- ALICE. YIELDX [44141] 


89.0 


4.8 


77.0 


14.0 


9.7 


3.0 


7.2 


5.2 


8.0 


Avignone [43] (MC) 


29.6 




34.4 


5.3 


4.4 




2.7 






Avignone 1431 (exp) 


30.0 ±7.0 




38.0 ±6.0 


3.5 ±0.9 


2.9 ±0.4 




3.3 ±0.8 






Klapdor-Kleingrothaus 1451 


58.4 


6.6 


79.0 


16.1 


10.2 




9.1 


1.6 


8.4 


Miley [46) 


26.5 


4.8 


30.0 


4.4 


0.5 










Barabanov |47l 


81.6 


2.9 
















Formaggio 1481 


29.5 






154.0 




5.4 


4.8 






ACTIVIA 


45.8 


2.8 


29.0 


8.5 


6.7 


2.0 


2.7 


1.6 


3.4 



Table 3 

Induced rates of isotope production Yj (kg -1 day -1 ) in enriched germanium (86% 7B Ge and 14% 74 Ge) at sea level. 





68 Ge 


60 Co 


65 Zn 


58 Co 


57 Co 


54 Mn 


63 Ni 


55 Fe 


HMS- ALICE. YIELDX |44l4l| 


13.0 


6.7 


24.0 


6.2 


2.3 


5.4 


6.0 


2.3 


Avignone 1431 


0.9 




6.4 


1.8 


1.0 


1.4 






Miley [151 


1.2 


3.5 


6.0 


1.6 


0.1 








Barabanov |47| 


5.8 


3.3 














ACTIVIA 


7.6 


2.1 


10.4 


5.5 


2.9 


2.2 


1.4 


1.6 



Table 4 

Induced rates of isotope production Yj (kg -1 day -1 ) in natural copper at sea level. For the ACTIVIA calculations, there are 
two sets of values; those from only the semi-empirical formulae and those that also include the MENDL-2P data. 





56 Co 


57 Co 


58 Co 


60 Co 


54 Mn 


59 Fe 


46 Sc 


IDEA [TT] 


22.9 


88.3 


159.6 


97.4 


32.5 


6.5 


3.8 


Baudis [49] 




30.5 




25.7 


134.2 






Formaggio |48j 


972.3 




131.5 


2675.5 


52.2 






ACTIVIA formulae only 


8.7 


32.5 


56.6 


26.3 


14.3 


4.2 


3.1 


ACTIVIA with MENDL-2P 


14.1 


36.4 


38.1 


9.7 


12.5 


1.8 


3.1 



and [4] are systematically higher than those from 
other analyses; it is unclear why this is the case. 

One important source of uncertainty in the pro- 
duction rates is the cosmic ray spectrum, where dis- 
crepancies around a factor of two can be explained 
by differences between the parameterisations used 
in various calculations (such as the spectrum shown 
in Fig. [2] used in ACTIVIA) and cosmic ray fluxes 
(which depend on geographic location and time) 
measured by various experiments. 

Another important issue is that the semi- 
empirical formulae must be able to predict the 
cross-sections of isotope production over a wide 
range of isotope masses and input beam energies. 
As shown in Fig. [TJ there are several physical pro- 
cesses that can take place in nuclear reactions. Each 
process must be described by a set of formulae that 
can adequately predict the cross-section for any 



given isotope and beam energy. Obviously, this can 
only be accomplished if there is enough experimen- 
tal data to derive such functions and obtain all 
relevant parameters with sufficient accuracy. The 
data used by Silberberg and Tsao did not cover 
all possible isotopes and energies, which limits the 
predictive power of the formulae for various cases. 
As can be seen in Figs. [4] and El some isotopes have 
very few measurements available. More experimen- 
tal data are needed so that improvements can be 
made to semi-empirical formulae, data tables and 
Monte Carlo simulations for isotope production 
calculations. 

7. Code modifications and extensions 

The code is structured so that modifications or 
extensions can be made to several components of the 
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cross-section and yield calculations whilst keeping 
the overall framework intact. 

Different input beam spectra can be used by 
including classes that inherit from ActBeamSpec- 
trum; the appropriate spectrum must be chosen 
by the input class. For example, a class describing 
the flux of an ion beam can be provided, although 
the formulae for calculating the cross-sections of 
products from such a beam would need to be im- 
plemented as well. The beam spectrum class also 
needs to specify the energy range and bin width 
that are used to calculate the integral in Eq. [21 as 
well as the atomic number Zb eam and mass number 
^4-bcam of the beam particles. 

The "decayData.dat" file contains a list of all 
available isotope products with half-lives greater 
than one day, as well as their side-branch nuclei. 
Other isotopes, with shorter half-lives, can be added 
to this data file. The code will calculate the cross- 
sections and yields for these nuclei automatically 
without any changes necessary to the source code. 
However, any side-branches matching the new iso- 
topes would need to be removed to avoid double 
counting. 

At present, the cross-sections are calculated by 
the ActSTXSecAlgorithm class using a combination 
of values from (MENDL) data tables as well as semi- 
empirical formulae from Silberberg and Tsao. Cross- 
section models for processes like spallation and fis- 
sion are implemented in their own classes, making it 
easier to update their parameters or formulae when 
necessary. New nuclear models can be added to the 
ActSTXSecAlgorithm class if required. 

Other cross-section algorithms can be specified, 
as long as they inherit from the base class ActAb- 
sXSecAlgorithm and are selected appropriately in 
the input class. For example, new algorithms may 
be needed to specify cross-section formulae when 
beams other than cosmic rays (protons/neutrons) 
are used. 

Other decay algorithms could be written, such 
as those that take into account chained decays, if 
they follow the Act AbsDecayAlgorithm abstract in- 
terface, which has access to the list of possible iso- 
tope products (ActProdNuclideList), as well as the 
definition of the target (ActTarget) and the Act- 
ProdXSecData results it stores (the target-product 
cross-sections). ActNuclide objects can be retrieved 
from the target and product list to take care of any 
bookkeeping required for chained decay algorithms. 
Note that the decay algorithm is separated from the 
cross-section calculations, allowing the possibility of 
using different decay algorithms for the same initial 
production yields, provided they are selected appro- 



priately by the input class. 

Changes can be made to how the calculations are 
set up. The default Actlnput class uses standard 
C ++ input/output commands to get the necessary 
information from the user. Other classes using dif- 
ferent input formats, such as XML files, could be im- 
plemented. However, they must all derive from the 
abstract ActAbsInput interface, and provide all re- 
quired inputs (beam spectrum, target and product 
isotopes, decay times and selection of cross-section 
and decay yield algorithms). 

At present there are two output formats available: 
ASCII text files and ROOT [IU\ files. Other output 
formats could be written, provided the classes follow 
the abstract Act AbsOutput interface and know how 
to deal with lines of text, tables of data (ActOut- 
put Tables) and graphs (ActAbsGraph objects). 

8. Summary 

We have presented the ACTIVIA C++ com- 
puter package that can be used to calculate target- 
product cross-sections and the production and de- 
cay yields of isotopes from cosmic ray activation. 
Cross-sections are calculated using a combination 
of semi-empirical formulae, assumed to be the same 
for proton and neutron beams, and ASCII tables of 
(experimental) data. 

The basic structure of the code, designed to work 
in the Unix environment, has been presented. It 
is structured so that modifications or extensions 
can be made to several components of the cross- 
section and yield calculations whilst keeping the 
overall framework intact. Since ACTIVIA uses 
only formulae or data tables to calculate cross- 
sections, and does not simulate individual particles 
and their interactions, it is easier to use than full 
Monte Carlo simulation computer packages such as 
GEANT4 [51] , while it is more extensible than sim- 
ilar FORTRAN-based codes like COSMO [42]. At 
present, the code is limited to only calculating the 
isotope production cross-sections and yields for one 
given target element (with its various isotopes); for 
targets made up of several elements, the code must 
be run separately for each and the results combined 
at the end by the user. Results from the calculations 
can be written out as ASCII text or stored within 
ROOT [TO] files. 

We have compared isotope production results 
from ACTIVIA against experimental data and es- 
timates from other calculations; in general, they 
agree within a factor of two. It is possible to change 
or extend the algorithms and formulae used to cal- 



culate the cross-sections and isotope yields. Finally, 
we have discussed ways of modifying the code such 
that it may be used in applications other than 
calculating yields of cosmogenic isotope products. 
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